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The interplay of deconfinement and chiral symmetry restoration are considered in terms of effective 
theories. We generalize the earlier model studies by considering fermions in higher representations, 
and study the finite temperature phase diagrams of SU(2) and SU(3) gauge theories with two 
fermion fiavors in fundamental, adjoint or two-index symmetric representations. We discuss our 
results in relation to recent lattice simulations on these theories and outline possible applications in 
the context of dynamical electroweak symmetry breaking. 

PACS numbers: 



I. INTRODUCTION 

Our understanding of low energy QCD, and properties of hot and dense matter and phase diagrams of strongly 
interacting theories in general, is rooted in exact and approximate symmetries. For QCD with light quarks, the 
spontaneous breaking of the chiral symmetry is manifested in the spectrum of hadronic states, while for pure gauge 
theory the properties of the phase transition between confinement and hot gluonic matter are captured by universality 
arguments based on the center symmetry. The perturbative analyses, although developed to impressive orders are 
inapplicable in the vicinity of the phase transition. Since this is the most interesting region for applications, like the 
heavy ion collisions, for hot QCD, alternative methods must be applied. The first principle method is provided by 
the lattice simulations and complementary methods are provided by the gauge/gravity dualities and effective models. 

During recent years, a simple model framework able to account for the observed behavior in QCD has been devel- 
oped. The model studies have mostly concentrated on hot and dense QCD, but also the dependence on the symmetry 
breaking parameters has been investigated. In this paper we will consider the Nambu-Jona-Lasinio model 2, i (NJL 
in the following, see [IHS] for reviews), enlarged with a coupling to the Polyakov loop [7HS] (PNJL model, see 
and references therein for previous studies), and extend it by considering higher representation matter fields. 

Generally, in a gauge theory, the matter fields in other than the adjoint representation break the center symmetry 
explicitly. In order to understand generic properties of strongly interacting gauge theories with matter in arbitrary 
representation, one needs to understand the fate of the deconfining phase transition as a function of fermion repre- 
sentations. For example, for SU(3) gauge theory with two Dirac flavors in the adjoint representation there are two 
transitions clearly separated: the deconfinement, measured by a rapid rise in the average value of the Polyakov loop, 
is observed at temperatures ca. an order of magnitude smaller than the onset of chiral symmetry restoration [26 . In 
comparison, for QCD with two flavors in the fundamental representation, only a single crossover with simultaneous 
rise of the Polyakov loop and decrease of the chiral condensate is observed. 

The main phenomcnological motivation for this study is in the possible applications of such theories in the model 
building for electroweak symmetry breaking via strong dynamics. As is well known, when coupled with the electroweak 
currents of the Standard Model, the chiral condensate of light QCD quarks breaks electroweak symmetry, and provides 
dynamically the masses for the electroweak gauge bosons. This scenario is, however, phenomenologically unrealistic: 
the masses of the W and Z bosons would be only ~ 30 GeV and, moreover, the triplet of pseudoscalar Goldstone boson 
appears in the physical spectrum rather than being absorbed into the longitudinal modes of the massive electroweak 
gauge bosons. Nevertheless, given the hierarchy and naturality problems of the elementary scalars, one can entertain 
the thought that in place of the scalar sector of the Standard Model, a dynamical mechanism of the type described 
above occurs in a new strongly interacting sector of the Standard Model; this is the vintage Technicolor introduced 
by Susskind and Weinberg in late '70s |271[55]. The original models, which simply replicate QCD-like dynamics, are 
mostly ruled out by the electroweak precision data |29| . Currently the model building efforts are concentrated on 
models which are quasi-conformal, meaning that the dynamics of the theory is governed by a quasi-stable infrared 
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fixed point. Over a large hierarchy of scales the Technicolor coupling constant, then, evolves very slowly, walks, in 
contrast to the running of the precociously asymptotically free QCD coupling. Hence these theories are categorized 
as walking Technicolor. The perturbative beta-function gives a guide how to construct such theories: the competition 
between anti-screening non-abelian gauge fields and screening matter fields implies a critical number of fiavors N* 
above which there exists a nontrivial infrared fixed point. Comparing the value of the fixed point implied by the two 
loop beta function and the critical coupling for the onset of the chiral symmetry breaking defines the lower boundary 
of the conformal window of the gauge theory under consideration. For walking Technicolor one wants to tune the 
matter content of the theory to reside close to the lower boundary of the conformal window, but still remain in the 
confining phase. 

For fundamental representation fermions the conformal window is expected to lie, roughly, above Nf -iNc- If 
all these fiavors are charged under the electroweak interactions in the usual manner, one typically encounters again 
large tension with the data from the electroweak precision measurements. This is so, since the so-called S'-parameter 
effectively counts the new degrees of freedom contributing to the electroweak symmetry breaking. A way to alleviate 
this tension was proposed in |30j : Since the fermions in higher representations carry more charge, they screen more 
efficiently, and as a consequence walking can be achieved with fewer flavors. 

In this paper we initiate a systematic study of higher representations within the PNJL model with the aim of 
establishing the phase structure, in particular the intertwinement of deconfinement and chiral symmetry restoration. 
We consider the strong dynamics in isolation, i.e. we do not consider the coupling with the electroweak currents. Also, 
to better illustrate the fixing of model parameters between theories with different fermion representations, we use as a 
benchmark model the PNJL model tuned to QCD, which will also allow us to compare our results with respect to the 
ones currently available in literature |31j . In this paper our aim is to show that theories with higher representation 
matter fields lead to novel patterns of how deconfinement and chiral symmetry restoration intertwine; we expect these 
patterns to remain intact if the theory is scaled to correspond with the strong dynamics at the electroweak scale. 
Quantitative applications for particular Technicolor models, which would also require the inclusion of electroweak 
currents, we leave for future work. 



II. MODEL FRAMEWORK AND PARAMETER SETUP 



A. Model setup 

In this paper we consider the Polyakov extended Nambu-Jona-Lasinio (PNJL) model. Based on our earlier experi- 
ence, we expect that similar results will be obtained also in the Polyakov extended quark meson models (PQM). The 
PNJL model is defined by the Lagrangian 

^ '^chiral "t" '^Polyakov '^interaction- (1) 

The chiral part of the Lagrangian corresponds to the usual NJL model. The pure gauge dynamics are encoded in 
'Cpoiyaiiov, which is essentially a parametric fit to lattice data. Finally, the interaction term features a connection 
between the chiral and pure gauge sectors. This is achieved by coupling the quarks to a static background gauge 
field; the derivation has been exposed in numerous earlier works, see e.g. [321133] . For the purposes of this paper it is 
sufficient to state the end result which is the grand potential of the form 

n = u^ + Ue + n^^. (2) 

The chiral part, U-^, is given by 

a, = '^^^ w 

where constituent quark mass is M = niQ — G{qq). 

The gauge sector, Ui, is independent of the chiral model but depends on the chosen gauge group. In this work we 
consider SU{2) and SU (3) gauge groups. The respective Z{2) and Z{2>) symmetric forms for the pure gauge potential 
are 

Ui,^,.,lT^ = -"f^KF,,<„|' + 6(T)ln[l-Kfs„(.,l'] (4) 
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where Ip the Polyakov loop in the fundamental representation of the gauge group, 
coefficients a(T) and b{T) are 



The temperature dependent 



a{T) = ao + ai ( ^ 



a2 



and b{T) = 63 



(6) 



where the parameters ag, ai, 02, and Tq will be fitted to pure gauge lattice data. The arguments of the logarithms are 
based on the Haar-measure of the gauge group [32] and results in a divergence at £p = 1, thus restricting the Polyakov 
loop to values £f < 1. In what follows, we will assume that the Polyakov loop expectation value is homogeneous 
in space. This assumption will simplify the evaluation of the fermion determinant; see below. The final part, f2qg, 
of the grand potential, contains the interaction terms between the chiral and gauge sectors. It arises from a simple 
integration over the quark fields coupled to a background gauge field 



^qq — 



^dim(i?)2iV/ 



(2^)3 



E 



-2NfT 



(27r)3 



Tr^ In 



1 + LjiC 



Tr„ In 



(7) 



where E — + and the constituent mass M is given below equation ([3|). Finite baryon chemical potential, 
/i, will not be considered in this paper, and in subsequent equations we assume that fi = 0. In the above equation, 
the dimension of the fermion representation, dim(i?), is made explicit. In the above equation the Wilson line in the 
representation R of the gauge group is defined as 



Lii{x) = Pexp 



l/T 



dTA^{x,T)T''{R) 



= exp[-^gr"(i?)/r], 

where the latter equality holds for a static background field. The Polyakov loop is then defined as the trace 

TrL_R(^)- 



dim(i?) 



(8) 



(9) 



The first term on the r.h.s. of Eq. (?]), which represents the vacuum fluctuations, is divergent and needs to 
be regulated. In PQM, the renormalizability of the model allows for a rigorous renormalization procedure of the 
divergence [M]. On the other hand, in the case of the PNJL model we will use a momentum dependent NJL coupling. 



G{p)=^G9{A~\p\) 



(10) 



In the vacuum term, this amounts to cut off the condensate dependent integral at the scale A, thus making the latter 
convergent (the divergence appears in a condensate independent integral which can be subtracted) . The thermal part 
of I?]) is convergent, hence it does not need regularization. However, the prescription ( [lO| implies that for p > A 
the free gas contribution M — rriQ has to be taken. These modes do not give an explicit contribution to the chiral 
condensate; however, they are Z3 charged and thus couple to the Polyakov loop. As a consequence, they contribute 
to the thermodynamics of the system and have to be taken into account. The prescription in ( 10 1 corresponds to an 
oversimplified version of a nonlocal and static NJL vertex. In our study, the choice (10) is justified a posteriori, since 



only in this case we obtain clearly separated deconfinement and chiral restoration in the case of adjoint fermions. 
This result is in agreement with what has been found in [31 where two colors QCD with adjoint fermions has been 
considered. As a matter of fact, in 31j a cutoff on the thermal part of the free energy density is introduced, in order 
to have a net separation of the two QCD transitions in agreement with the lattice results. 

The thermodynamics of the system is obtained by solving the equations of motion for the order parameters of the 
system 



on 

dM 



-0, 



W) 



= 0, 



(11) 



and then evaluating the grand potential at the minimum to obtain the free energy at given temperature. We work 
under the mean field approximation, replacing the Polyakov loops with their thermal expectation values, £ — >■ (£). 
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B. Introducing higher representation fermions 



Our main objective is the study of fermions in higher representations. Besides these, we consider fermions in the 
fundamental representations of the color group in order to fix the model parameters. For the SU(2) gauge group we 
use the fundamental and adjoint representations while for the standard SU(3) case we include also the two-index 
symmetric representation, i.e. the sextet, in our study. As explained in the Introduction, these cases are relevant for 
the phenomenology of walking technicolor theories. In this section we will describe how we construct the fundamental 
and higher Polyakov loops and give the form of the interaction potential (7| for the different representations. 

In the fundamental representation of a SU{N) gauge group, the Wilson line Lp can be written as a diagonal matrix 

L;^=diag(e*''\--- ,e^^"), (12) 

with the constraint 9n = + • • • + On^i). For SU{2) this means that there is only one free parameter 9; for the 
SU{3) group there are two independent parameters, 9i and 02- For simphcity, in the SU{3) case we choose ^2 = 
which leaves us with one free parameter 9i = 9. This assumption is justified at zero quark chemical potential, which 
is the case we consider in this article, since the Polyakov loop is real, £ — £'''. 

The Wilson line for the adjoint representation (for both SU (2) and SU{3)) can be constructed from the fundamental 
Wilson lines using the relation 

Lf = 2Ti[LfT''L^pT''], (13) 

where r"s are the generators normahzcd as Tr[rT''] S''''/2. In the SU{3) case, in addition to the adjoint Wil- 
son line, we want to know the Wilson line in the two index symmetric representation. The Wilson line for this 
representation is related to the fundamental ones through a similar relation 

Lf = TiiLFT^LFT^], (14) 

where the set of matrices matrices Tg (i = 1 • • • 6) is a basis for symmetric 3 x 3~matrices and given explicitly in 
Appendix 

The Polyakov loops can be calculated as a function of 9 according to the definition ([9|. For the fundamental 
representation one gets 



%.(3, - :;(l + 2cos0). (16) 



^Fsu(2) = cos 61 (15) 

1 

3* 

The higher representation loops can be similarly calculated and expressed in terms of the fundamental loops 

= ^(2cos20 + l) = i(44^^^^^-l) (17) 

^Asu.s, - ^(2 + 4cos^? + 2cos20) = ^(94^.^^3^-l) (18) 

4.,(3) = ^(2 + 2cos^? + 2cos20) = ^(94^.^^3^-3^^,,,3,). (19) 

One can now proceed to perform the color traces in potential ([t]) for the different gauge groups and representations. 
For this we write the interaction part of the potential as 

n,, = -ANfT J ^ (Sj,) dim{R)2Nf J ^E, (20) 



where we have defined 

5fl, = Tr,ln[l + L^e-^/^]. (21) 

Performing the traces in SU(2) we find 

Sf,^,,^ = ln[l + 2£f,^^^^ e-^l^ + e^^i^/Tj (22) 

and 

= ln[l + e-^l^\ + ln[l + (4£^,„,,, - 2)e-^/^ + e'^^l^]. (23) 
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TABLE I: Parameters of the NJL model and derived observable values for different cases. The SU(3) fundamental values 
correspond to observed QCD properties, while the rest are obtained by scaling the coupling G as explained in the text. 





NJL parameters 


Derived observable values 




mo G A 


m,r nicr M 


SU(3) Fundamental 
SU(2) Fundamental 
SU(3) Adjoint 
SU(2) Adjoint 
SU(3) Sextet 


5 MeV 8.879 GeV"'' 678 MeV 
5 MeV 13.32 GeV'^ 678 MeV 
5 MeV 7.492 GeV^ 678 MeV 
5 MeV 23.66 GeV'^ 678 MeV 
5 MeV 11.099 GeV-^ 678 MeV 


93 MeV 138 MeV 599 MeV 296 MeV 
76 MeV 137 MeV 598 MeV 296 MeV 
167 MeV 166 MeV 2567 MeV 1157 MeV 
96 MeV 179 MeV 2811 MeV 1406 MeV 
139 MeV 174 MeV 2615 MeV 1304 MeV 



TABLE 11: Lattice data fitted coefficients for the SU(3) and SU(2) pure gauge potentials ( [5| a nd Q, respectively. Values for 
SU(3) are obtained from 9^, while SU(2) values are fitted in this worlt to lattice data from [37]. 



SU(3) 


ao 


ai 




bs To 


3.51 


-2.47 


15.2 


-1.75 270 MeV 


SU(2) 


ao 


ai 


a2 


bs To 


1.19 


-1.136 


7.94 


-2.759 268 MeV 



Analogously, in the case of SU(3) we find 

^;^,„,3, - ln[l + e-^/^] + ln[l + (Mp,,,,^ - l)e-^/^ + fi-^^/^], (24) 



5a,„(3, = 2 ln[l + e-^/^] + 2 ln[l + (3£j.,,,3, - l)^-^/^ + g-^^/^] 

+ ln[l + (94,,,3) - 6%.(3) - + e-'^n ^ ^ 

and finally 

56.„,3, = 21n[l + e-^/^] + ln[l + (O^^^^^^ - 3^^,„,3, - 2)(e-^/^ + e~'^/^) 

We have written the color traces for the different representations in terms of the fundamental representation loop 
of the corresponding gauge group. Proceeding this way we keep the same observables, the constituent mass and 
the fundamental Polyakov loop, when comparing the different representations and their effects on the chiral and 
deconfinement transitions. 



C. Fixing the parameters 

In this section we discuss how the model parameters, fixed for three colors and fundamental representation fermions 
can be used to obtain the parameters in the case of higher representation matter or two colors. There are three 
parameters which enter into the chiral sector of the PNJL model. These are the bare quark mass mo, the value of the 
four-fermion coupling G and the momentum cut-off A. For the SU(3) fundamental fermions, we fix the parameters 
as usual [H |51 1351 136j by requiring that the physical values of TOtt , /^r and (qq) are reproduced. The parameter set we 
use is shown in the first row of Table HI 

In addition to the parameters of the chiral sector, the parameters for the pure gauge potentials ([5| have to be 
assigned. This can be done by fitting the potential to pure gauge lattice data. In this work we will use for SU(3) the 
fit presented in [9 and the required coefficients are shown in Table [ll| 

With the set of parameters in Tables [l] and [llj the thermodynamic potential for the PNJL model in the case of 
SU(3) with fundamental quarks is completely specified. 

To discuss higher representation fermions and different numbers of colors, we need a scheme to translate the model 
parameters between different cases. As a matter of fact, in the case of representations different from the fundamental 
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TABLE III: Quadratic Casimir operators and representation dimensions for SU(N). 



Representation 


Ca 


dim 


Fundamental 


JV^-l 
2JV 


TV 


Adjoint 


N 


iV^ - 1 


Symmetric 


[N-l)(N+2) 
N 


JV(JV+1) 
9. 



one of SU(3), the numerical values of physical observables are not known. As a consequence, we need some argument 
to obtain the parameters for the higher representations, from those fixed for the fundamental representation. We 
will first discuss how the parameters are obtained for higher representations of SU(3) and then how to obtain the 
parameters for the two color theory. 

Of the three parameters, toqi G and A, the quark mass is basically a parameter that determines the amount of 
explicit chiral symmetry breaking in the theory and is thus responsible also for the mass of the pion. It is essentially 
a free parameter even in QCD so there is no need to adjust it when changing fermion representations. The remaining 
parameters G and A are connected to the fundamental theory, in the case under consideration to QCD with two colors 
and/or fermions in different representations. For what concerns the cutoff, we do not change its numerical value after 
changing the fermion representation. A more rigorous treatment would require the fixing of A (and of G as well) by 
reproducing the phenomenological properties of the vacuum as in the case of three colors discussed above. However, 
data about these properties are missing. In absence of a serious guiding principle, we use the same value of A for all 
the cases considered in this article. 

On the other hand, we have a theoretical guiding principle to change the value of G according to the dimension of 
the representation considered, which was introduced in |31j . It is based on the Fierz transformation properties of the 
color current interaction 



A„t = -.gs(V57^T;flV)' (27) 

in which the coupling gs can be related to the underlying QCD coupling and the NJL coupling G For an effective 
NJL interaction of fermions in representation R, 

-CNJLint = -Gniqq)^ , (28) 



one gets 

where G2{R) is the quadratic Casimir operator for representation R. For the derivation of this result we refer to the 
original article [3T]. Assuming that the coupling gg is the same for all representations, the ratio of the NJL couplings 
for different representations will be independent of gs '■ 

Gr, ^ C2{Ri)dimR2 

Gr, C2iR2)dimRi ■ ^ ' 

By means of the above equation we are able to determine the higher representation NJL coupling from the fun- 
damental one, once the Casimir operators for the relevant representations of SU(N) are specified. These operators 
for the fundamental, adjoint and symmetric representations are shown in Table |III[ Together with the fundamental 
representation ones, the model parameters for higher fermion representations are collected In Table |l] 

Then, let us consider setting up the model parameters in SU(2) case. As in the SU(3) case with higher representation 
matter, we have no experimental data to fit that would determine our model parameters. Therefore we again make 
use of scaling arguments to infer the SU(2) parameters from the known SU(3) ones. In particular, in the large Nc 

1/2 

limit of QCD (see [38] for a review) one finds f-^ Nc and (qq) cx Nc- These scalings are respected by the NJL 
model, see for example [1]. Therefore it is natural to fix the values of G and A, for fermions in the fundamental 
representation of SU(2), by requiring that the values of the properly scaled and {qq) are reproduced. For what 
concerns the value of the bare quark mass, the GMOR relation f^m^ oc ■mo{qq) together with the scaling cx 
implies tuq oc N^, hence mo does not need to be changed. 

In the NJL setup used here, this is consistent with assuming that the four-fermion coupling G is proportional to 
l/Nc- To see this, one may first note that the NJL gap equation, with cutoff A ^ N^, implies M ~ N^. Then, since 
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the gap equation can be written as AI = mo + 2G{qq), one further obtains {qq) ~ Nc consistently with the large Nc 
expectation discussed above. For the determination of /^r and m^r are determined via Eqs. (4.26) and (4.21) of [1], 
respectively, and these imply that /^r ^ Nc and ~ again consistently with the expected large Nc behavior. 

Thus going from fundamental SU(3) to fundamental SU(2) we multiply G by 3/2 and, to keep things as simple as 
possible, we again do not alter the momentum cutoff. As a result of this scaling choice the pion and sigma vacuum 
masses, calculated from the NJL model, remain the same between fundamental SU(3) and fundamental SU(2). Then, 
to define the model corresponding to SU(2) with two adjoint fcrmions we use the same Fierz scaling argument that 
we utilized in the SU(3) case. The resulting parameter values for both fundamental and adjoint SU(2) are also shown 
in Table U 

Finally, since also the pure gauge sector is changed when switching from SU(3) to SU(2), new parameters have to 
be assigned for the pure gauge potential Q as well. As in the SU(3) case the coefhcients arc fitted to reproduce 
lattice data. We use SU(2) pure gauge results from [35] and fit the coefficients of the potential in Eq. Q so that it 
reproduces a matching energy density as a function of temperature with the data. The data with the resulting fit are 
shown in Figure [l] while the numerical values of the coefficients ao, ai, 02, and Tq are listed in Table [TTl 




2000 2500 



T [MeV] 



FIG. 1: The data points show the energy density of SU(2) pure gauge theory [37] > while the solid curve shows the corresponding 
result obtained from the Polyakov loop potential Q with parameters as given in Table In] 



III. NUMERICAL RESULTS 

In this section we establish the finite temperature phase diagrams of PNJL model with fermions in fundamental and 
higher representations of the SU(2) and SU(3) gauge groups by numerically evaluating the temperature dependence of 
the Polyakov loop condensate and the constituent mass. These are the (approximate) order parameters for the chiral 
and deconfinement transitions, and their temperature dependence displays the interplay between the deconfinement 
and the chiral restoration crossovers. 

First, we report on our results on the SU(2) gauge group. This case has been already studied in the literature 
by means of effective models coupled to the Polyakov loop [31] [40] ; therefore the results presented here have to be 
considered as a check of our calculations. Nevertheless, they are quite important since they are useful to interprete the 
analogous results obtained in the case of the SU(3) gauge group. Before going ahead, we notice that in the numerical 
computations with fermions, we do not change the value of Tq in the Polyakov loop effective potential. This is done 
just for a matter of simplicity; in principle, an explicit dependence of Tq on the number of active flavors should be 
introduced [41] . This choice for simplicity, however, allows us to focus only on the few parameters in the chiral sector, 
which will be changed according to the representation and to the gauge group as explained in the previous subsections 
(at the end of the day, we just will modify the coupling constant). This is helpful to realize what is the main source of 
the different behavior of deconfinement and chiral symmetry restoration when different representations for fermions 
are considered. 

Finally, we ignore the explicit dependence of the NJL coupling on the Polyakov loop, as studied in the recent 
literature [2TJ l42l - l44] . Also this choice is inspired by simplicity. Together with other compromises we have made 
in fixing the values of the model parameters, also this aspect has a quantitative impact on the interplay between 



8 



deconfinement and chiral symmetry restoration, and certainly it deserves further study which we leave for a future 
project. 

For the SU(2) case, the chiral and deconfinement order parameters as a function of temperature are shown in 
Figures [2] and [3] for the fundamental and adjoint fermions, respectively. In the case of fundamental fermions we expect 
the behavior of the order parameters be similar with the more familiar SU(3) fundamental fermion case; indeed both 
chiral and deconfinement transitions are smooth cross-overs which are located quite close together. However, for a 
crossover transition the exact transition temperature depends largely on the definition one uses. From Figure [2] one 
could say that both transitions happen near T ~ 220 MeV since near this temperature both the normalised constituent 
mass and the Polyakov loop cross the value 0.5. If one makes a comparison with the respective SU(3) case, shown in 
Figure |4j one can note that perhaps the only significant difference is the slightly reduced transition temperature and 
a smoother rise in the value of the Polyakov loop when approaching the phase transition region T = 200 — 220 MeV. 
This close similarity can also be expected by looking at Table [T] where the pion, sigma and constituent quark masses 
show no deviation between the SU(2) fundamental and SU(3) fundamental cases. As explained earlier, this is due to 
the scaling of G being partially cancelled by the change in the number of color degrees of freedom. 

In the case with adjoint fermions the situation is different from the fundamental case. A large separation between 
the two transitions appears, and the deconfinement transition has changed qualitatively. This is in agreement with 
earlier studies [iljl and with the lattice results. In more detail, the pseudocritical temperature of the chiral cross-over 
is much higher than the one found in the case of fundamental fermions, namely around 700 MeV. On the other hand, 
the deconfinement transition is now a true second order phase transition near T 175 MeV. The latter aspect is 
easily understood: Fundamental fermions break the center symmetry, thus act as a source for the Polyakov loop, and 
Z(2) breaking persists in the confining phase, turning the deconfinement transition to a crossover. On the other hand, 
adjoint fermions do not break Z(2), which turns to be a true symmetry of the action, unbroken in the confinement 
phase. This is what one would expect solely on the basis of underlying symmetries |45j . 




800 



T [MeV] 



FIG. 2: The chiral and deconfinement order parameters as a function of temperature for the PNJL model in the SU(2) 
fundamental fermion case. Model parameters listed in Table |l] 

Let us then turn to the numerical results on SU(3) theory with fermions in fundamental, adjoint or sextet repre- 
sentation. As in the SU(2) case, we consider the behavior of the order parameters as a function of temperature. The 
results are shown in Figures |4j [5] and [6] for the fundamental, adjoint and sextet representations, respectively. 

The order parameters for the case of fundamental SU(3) fermions show only minor quantitative differences with 
the respective SU(2) case and feature the typical coincidence of the two crossovers in a range of temperature centered 
on Tc « 230 MeV. Note that this crossover temperature is a little higher than in the earlier works |331llSlll7j. The 
difference is small, ^ 5 percent, and is not of importance for the results of this work. The discrepancy is due to 



two reasons. First, the dynamical mass cutoff in thermal contribution of quarks, arising from Eq. (10), was not 
accommodated in the calculations of [331 SSI EZ] • Second, we are using here slightly different initial parameter set 
than in [531 HSl HZ] , but its effect on the transition temperatures is very small. 

As in the SU(2) case, the temperature dependence of the order parameters for adjoint fermions of SU(3) correspond 
to expectations based on the approximate symmetries of the order parameters. The NJL model exhibits separate 
transitions with roughly a 400 MeV gap between the two. Again the chiral transition is a smooth crossover while the 
Polyakov loop retains the first order transition from the pure gauge sector. The deconfinement, which is a crossover 



9 




FIG. 3: The chiral and deconfinement order parameters as a function of temperature for the PNJL model in the SU(2) adjoint 
fermion case. Model parameters listed in Table |l] 



for fundamental representation matter fields, is turned to a true phase transition because quarks in the adjoint 
representation do not break the center symmetry Z(3) explicitly. 

Finally, we turn to the more interesting case of the PNJL model with sextet fermions, which has not yet been studied 
in the chiral model literature. The SU(3) sextet results are remarkably similar with the SU(3) adjoint case: both phase 
transition temperatures deviate less than 20 percent from the adjoint, shifted towards higher temperatures, and the 
qualitative features of the phase transition are nearly identical. At a first sight, this behavior of the order parameters 
might appear unexpected, since the center symmetry is explicitly broken by fermions in the sextet representation. 

Our interpretation of the results is as follows: In the case of sextet fermions, the constituent quark mass at zero 
temperature is considerably higher than the mass of fermions in the fundamental representation. This is the natural 
consequence of having a larger NJL coupling in the sextet representation. When quarks couple to the Polyakov loop, 
the amount of breaking of the center symmetry is proportional to the hopping parameter, which in turn scales as 
1/m where m is the quark mass. In the model at hand, at the one-loop level the mass eigenstates which propagate 
in the Polyakov loop background have mass m = M with M being the constituent quark mass. The latter is very 
large in the case of the sextet fermions, which causes the hopping parameter to be very small in this case. Thus, the 
explicit breaking of the center symmetry, even if present, is very soft, and the Polyakov loop behaves effectively as if 
the center symmetry was unbroken. 




T [MeV] 



FIG. 4: The chiral and deconfinement order parameters as a function of temperature for the PNJL model in the SU(3) 
fundamental fermion case. Model parameters listed in Table |l] 
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FIG. 5: The chiral and deconfinement order parameters as a function of temperature for the PNJL model in the SU(3) adjoint 
fermion case. Model parameters listed in Table |l] 




FIG. 6: The chiral and deconfinement order parameters as a function of temperature for the PNJL model in the SU(3) sextet 
fermion case. Model parameters listed in Table |l] 



IV. CONCLUSIONS 



In this paper we have initiated a study of PNJL model for apphcations beyond QCD. In particular we concentrated 
on the SU(2) and SU(3) gauge theories with two Dirac fermion flavors in higher representations of the gauge group, 
as these have been established as a phenomenologically viable theories for dynamical electroweak symmetry breaking. 
Eventhough these theories are interesting for understanding general features of strong dynamics, the real motivation 
for our study is that if such strong dynamics underlies the Higgs sector to be uncovered at the LHC, it becomes 
immediately of interest to study their finite temperature phase diagrams for the evolution of the early universe. 
Based on the experience gained in investigations of hot and dense QCD matter, the PNJL model provides a good 
quantitative tool for such studies. As a novel new result we have shown that, while the theories with fundamental 
and adjoint fermions show the behavior expected by considering the underlying symmetries, the SU(3) theory with 
two sextet fermions is, despite the explicit breaking of the center symmetry, akin to the theory with adjoint fermions 
where the center symmetry is not broken by the presence of matter fields. 

One can imagine several further avenues for research. First, all the theories with higher representation matter 
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we have considered are currently investigated in lattice simulations, see e.g. |48H51j . Similarly as with the case of 
QCD, for which the PNJL model has been successfully applied, for these theories we expect the future developments 
and increased precision in the lattice analyses of these theories to help in better fixing the model parameters. For 
applications to electroweak phase transition in the early universe one should couple the electroweak currents. For 
theories with only chiral fields see [S2] for recent developments. We aim to proceed towards these directions taking 
also into account the nontrivial behaviors due to the deconfincmcnt transition which we derived in this paper. 
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Appendix A: A basis for symmetric 3x3 matrices 

To define the Polyakov loop in the sextet representation, we use the following basis for symmetric 3x3 matrices: 

1/^ 

Tgi = ( I 7| = ( 1/^/2 




Tg^ = I 1 I T^ = \ 1/V2 

1/V2 

1/V2' 
T| = I I r| = I 




1/V2 
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